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ABSTRACT 

In this study we predict the total distributions of powerful (FR II) active double-lobed radio 
galaxies and ghost sources, and their observable distribution in the X-ray sky. We develop an 
analytic model for the evolution of the lobe emission at radio and X-ray energies. During jet 
activity, a double radio source emits synchrotron radiation in the radio and X-ray emission 
due to inverse-Compton (IC) upscattering by 7 ~ 10 3 electrons of the cosmic microwave 
background. After the jets switch off, the radio luminosity (due to higher 7 electrons) falls 
faster than the X-ray luminosity and for some time the source appears as an IC ghost of a 
radio galaxy before becoming completely undetectable in the X-ray. With our model, for one 
set of typical parameters, we predict radio lobes occupy a volume fraction of the universe 
of 0.01, 0.03, 0.3 at z = 2 (during the quasar era) of the filamentary structures in which 
they are situated, for typical jet lifetimes 5 x 10 7 yr, 10 8 yr, 5 x 10 8 yr; however since the 
inferred abundance of sources depends on how quickly they fall below the radio flux limit 
the volume filling factor is found to be a strong function of radio galaxy properties such as 
energy index and minimum 7 factor of injected particles, the latter not well constrained by 
observations. We test the predicted number density of sources against the Chandra X-ray Deep 
Field North survey and also find the contribution to the unresolved cosmic X-ray background 
by the lobes of radio galaxies. 10-30 per cent of observable double-lobed structures in the 
X-ray are predicted to be IC ghosts. The derived X-ray luminosity function of our synthetic 
population shows that double-lobed sources have higher space densities than X-ray clusters at 
redshifts z ^ 2 and X-ray luminosities above 10 44 erg s _1 . 
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1 INTRODUCTION 

We investigate how much of intergalactic space is filled with cur- 
rent and old radio lobes. Most gal axies seem to ha ve a supermassive 
black hole at their centres ( Mag orrian et al] 19981) , which can serve 
as a central engine for jet ejection. It is possible that most massive 
galaxies have at least one outburst of jet activity leading to a giant 
radio source in their lifetime (lasting ~ 10 8 yr), probably between 
redsh ifts 1.5 to 3 during the quasar era l lGopal-Krish na & Wiita 
l200lh . The radio galaxies will have extended emission in the ra- 
dio as well as extended X-ray emission in the keV energies due to 
inverse-Compton (IC) upscattering by 7 ~ 10 3 electrons of pho- 
tons that comprise the cosmic microwave background (CMB). The 
IC and synchrotron losses will downshift the higher energy elec- 
trons (7 ~ 10 4 ) responsible for GHz synchrotron radiation more 
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quickly than the lower-energy electrons that give rise to the X-ray 
emission. Thus for some period of time after the jet is switched 
off the source will appear as an inverse-Compton ghost of a ra- 
dio source before becoming completely radio and X-ray dark. The 
extended X-ray source HDF 130 (z — 1.99) has been interpreted 
to be such a radio galaxy with its jets tur ned off and only show- 
ing a double-lobed structure in the X-ray iFabian et al1l2009h , and 
the models developed in this paper are applied to that object in 
a companion paper. The CMB energy density is proportional to 
(1 + z) 4 , cancelling the dimming due to distance, and so ex- 
tended X-ray emission may be observable at both low and high 
redshifts. Previous work on the extended X-ray emiss ion of radio 
galaxies range from as early as iFelten & R ees (1969), which de- 
scribed the cosmic X-ray background (CXB) being attributed to 
a Compton-blackbody proc ess acting in double radio sources, to 
iGopal -Krishna et "all Jl989h . which is an analytic treatment of the 
time evolution of the radio output of ex panding radio lobes incor- 
porating IC losses, to ICelot ti & FabiarJ d2004l) . which considered 
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the relative abundance of double radio galaxies to X-ray clusters, 
and recently to the work bv lNathl JioToh which developed a simple 
model for the evolution of the X-ray luminosities of double radio 
sources. 

We create an analytic model accounting for injection of rela- 
tivistic particles into the lobes of radio galaxies and adiabatic, syn- 
chrotron and IC energy losses to describe the evolution of emission 
in the radio and the X-ray bands of a double-lobed radio galaxy 
which turns off after a typical jet lifetime. We can describe the to- 
tal population of active (i.e. jets switched on) and no longer active 
(jets off) sources from an empirically-inferred birth function and 
jet energy distribution. Of interest are the volume filling factor of 
the universe or filaments by radio lobes, their observable number 
density in X-rays, their contribution to the unresolved CXB and the 
X-ray luminosity function (XLF) of these extended X-ray sources. 
We explore and find that individual radio galaxy properties such 
as the index and minimum energy of the injection spectrum, and 
the ambient density, can significantly affect the time at which ra- 
dio sources fall below a given flux limit and consequently have a 
significant impact on the actual observed distribution of sources. 

We consider the distribution and evolution of pow erful Fa- 
naroff and Riley Class II (FR II) jFanaroff & Riievll 19741) sources 
to understand what fraction of filaments may be filled with radio 
lobes. We consider a cosmology with Ho — 70 km s _1 Mpc" 1 , 
fi M = 0.27, = 0.73. 



2 THE MODEL 

In order to study the evolution of the lobes of FR II radio galax- 
ies, we create an analytic model where the lobes are fed with syn- 
chrotron emitting electrons from a compact hotspot which undergo 
adiabatic expansion and synchrotron and IC losse s . We follow 
the fo r malisms found in iKaiser & AlexandeJ i 19971) . I Kaise r et al.l 
dl997l) . lBlundell et alj Jl999h and lNathlfeOlOh . We are interested in 
studying the evolution of X-ray and radio emission of these sources 
to beyond the jet lifetime tj . 

We pick one set of typical model parameters (case [A]) de- 
scribing the ambient environment and injection spectrum based o n 
parameters used in previous works bvlKaiser & AlexandeiU 19971) . 
IKaiser et al.l Jl997l) . lBrundell et al.l Jl999h and lNathl d20ld) . How- 
ever, we later show the importances individual parameters may 
have in altering the inferred distribution of sources (§[7}. Case [A] 
consists of parameters inferred from observations and similar to 
parameters used in previous models for radio sources, and so it 
serves as a normative case. However, not all parameters are well 
constrained by observations, such as j m i n (see discussion on this 
in § 12.61 1, and so some of the other cases from [B]-[L] may also be 
quite plausible models for a typical radio lobe. 



2.1 Synchrotron radiation basics 

The synchrotron luminosity <5P„ in units of W Hz -1 sr -1 , aver- 
aged over pitch angle, of a volume element SV with magnetic en- 
ergy density ub is given by: 



1 Y 
6P V = — ct t cmb — n0j)5V 

07T V 



(1) 



where j is the Lorentz factor associated with an electron radiating 
at frequency v, 71(7)0(7 is the number density of electrons with 
Lorentz factors between 7 and 7+d7, ctt is the Thomson scattering 
cross section and c is the speed of light. We assume that electrons 



only emit at their peak frequency given by v — 7 v^, where is 
the Larmor frequency. 

The energy loss equation describes the time evolution of the 
Lorentz factor: 



(hi 



= -7 



1 1 dVi 



(7t 



3 V\ dt 3 m e 



7 2 (UB + Mc) 



(2) 



where the first term is the energy loss due to the adiabatic expan- 
sion of the lobe of volume V\ and the second term describes the 
synchrotron and IC losses. Here m c is the mass of an electron and 
2.7 K (1 + z)f is the CMB photon energy density at 
the redshift of the source, a — 7.565 x 10 -16 J K -4 m -3 is the 
radiation constant. 

Supposing we know the distribution of electrons 71(71, ti) as 
initially injected into the lobe at time ii and assuming the electrons 
are uniformly distributed over volume, we can find the distribution 
of those electrons at a later time, 71(7, i), from: 

n( 7 ,t)d 7 = n( 7 i,ti)^^d7 0) 
where determining ji and d r yi/d'y is described in § 12.31 



2.2 Dynamics of the radio lobes 

The density of the environment surrounding a radio galaxy will 
affect both the dynamics and energy evolution of the jets and lobes. 
We assume a density profile around the AGN described by 



p(r) = po(r/a ) 13 



(4) 



where for case [ A] we take = 1.5, a p = 10 kpc and po — 1-67 x 
10~ kg m " 3 i lBlundelletal]|l999l) . inferred from observations, 
but will also consider steeper and less dense profiles (§0. It may 
be the case that the radio-source environments change with redshift, 
but we will assume independence from 2 as the evolution of radio 
source environments is not yet fully understood. 

We assume in our model a double jet system with origin at 
the central engine of an AGN, each jet transporting a total power 
Qj, assumed to be constant while the jet is on. The jet will have 
a lifetime of tj , after which the jet activity stops and no further 
particles are injected into the lobe. During a time t < tj, we 
use the source expansion based on the characte ristic length scale 
(£ 3 Q j p " 1 a " ,3 ) 1/(5 - /3) described in lFalld Jl99ll) : 



Cl 



1/(5-/3) 



(5) 



where Lj (t) is the len gth of a jet and ci is a dimensionless constant 
taken to be 1.8 as in lBlundell et alj jl999h . The total length of the 
source is 2Lj . 

We make the assumptions o f lBlundelletal] dl999h that the jet 
remains relativistic from the central engine to the jet shock (i.e. 
hotspot), that the jet has lower density than the surrounding en- 
vironment, that at the outermost edge of the source at Lj(t) is 
the hotspot just beyond the shock structure where bulk kinetic en- 
ergy from the jet is randomized feeding the larger head (outermost 
bright emission region) of the lobe. The head pressure is responsi- 
ble for the rate at which the source grows. The hotspot pressure is 
responsible fo r adiab atic expansion losses out of the hotspot (see 
iBlundell et ail dl999h ). 

While the jet is on, the pressure in the head can b e found from 
the jump conditions of a strong shock, as discussed in Kaiser et al.1 
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we have: 



Phead(Z) 



2-/3 



(r x + i)(5-/3)- 



2 

:"•() 



Qi 



4+/3 



(6) 



where T x = 5/3 is the adiabatic index of the surrounding IGM. 

The lobe is treated as consisting of small volume elements S V, 
each with varying fluid properties. A fluid element is injected at 
some time ti from the hotspot over a time interval Sti. The element 
SV(U) is related to Sti by: 



SV{ti) 



(ri - i)Qj 



pi(ti) \cA hs pi(ti) 



6U 



(7) 



using thermodynamic relations assuming adiabatic expansion of 
the volume e l ement over the time interval (equation (19) in 
Blun dell et all (l999)). ylhs is the area of the hotspot (the hotspot 
is assumed to have a fixed radius of 2.5 kpc), T\ = 4/3 is the adi- 
abatic index of the lobe and p\ is the pressure in the lobe, which we 
equate with the pressure in the head given by equatio n ((6} divided 
by a factor of 6, adopted from Blun dell et alj dl999r) , to allow for 
the pressure gradient along the lobe. Then, assuming further adia- 
batic expansion, the volume element changes as: 



5V(t) = 



(Ti - l)Qj 
Pi(ii) 



cA hs pi{tC, 



»(*) 



Sti. (8) 



The volume of each lobe at time t can be found by integrating 
equation l[8j. That is, for t ^ tj : 



Vx(t) 



* (r t - i)Qj 



cA hs pi(ti) 



Pi ftp 
»(*) 



(9) 

For a general time t, we want to set the integration limit to 
min[t, tj], as nothing is injected into the lobes once the jet has shut 
down. 

We call the axial ratio of the lobe 7? — R(t) (nucleus-hotspot 
distance divided by full width of the lobe), so that the volume of a 
lobe is: 



Vi(t) = 



4R(ty 



-Lj(t) 3 



(10) 



(the lobe is assumed to have a cylindrical shape). For t ^ tj, we 
know how Vi and Lj evolve and thus can determine how R grows 
with expansion. 

Now we will describe the evolution of the lobe after the jet has 
turned off. The expansion of the lobe for t > tj is governed by 



(11) 



where p a is the ambient density at Lj(t), which we know from 
equation 10 • 

It is not reasonable to assume that the source continues to ex- 
pand according to equation ([5} once the jet activity has discontin- 
ued since at this point Qj = and we do not have a characteris- 
tic length. To proceed, we assume that the axial ratio for t > tj 
will be given by its value at R(tj), that is, the axial ratio at the 
time when the jet turned off. This is a reasonable assumption if the 
axial ratio is small enough and the distance from the galaxy large 
enough so that the external pressure at the heads and sides of the jet 
are similar, even in an ambient environment with density described 
by a power-law decline. The assumption is useful for obtaining an 
analytic solution for the problem, and has been made for the en- 
tirety of the evolution of sources in other models such as those of 



iKaiser & AlexandeJ dl997h . lKaiser et alj l fl997l) . and lNathl fcOld) . 
Noticing from equation ([9} that for t > tj we have 



Vi{t) = V x {tj) 



Pi(*j) 
Pi(*) 



(12) 



and using equation dlOt , we find after substitution from equa- 
tion Gj) that for t > tj , 



Lj = k^{L] 



•3IV2+/3/2x 



where 



2\r, 



(13) 



(14) 



With the assumption that R(t) = R(tj) for post-jet conditions, we 
find an analytic solution for the evolution of Lj at t > tj by solving 
the differential equation {T3} with appropriate boundary condition 
at time tj : 



h{t) = ((3Ti/2 + 1 - P/2)(Vk(t -tj) + O) 



2/(3^+2-/3) 



with 



C : 



Lj(tjf F ^ 2+1 -^ 2 
3Ti/2 + 1-/3/2 



(15) 



(16) 



From Lj(t) we can determine Vi(t) and consequently p\(i) 
from equation dl2b . for t > tj. 



2.3 Radio luminosity of the lobe 

We assume that the initial electron energy distribution when in- 
jected into the lobe is a power law in energy given by: 

n(7i,ti)d7i = n ^ p dy l (17) 

with, for case [A], p — 2.14 as in Kai ser et"al] J 19971) and 71 
ranging between 7mm = 1 and 7 max = 10 6 (Lorentz factors of 
7 ~ 10 3 are required to produce upscattering of the CMB in the 
X-ray and Lorentz factors of 7 ^ 10 4 are needed for GHz syn- 
chrotron radiation in the radio for typical magnetic field strengths). 
The constant no is found by integration: 



n 



U e (ti) 



(75 - l)7i P dji 



(18) 



Increasing 7 max h as little effect on the lumin osity. However, p can 
range from 2 to 3 (Alexande r & Leahvl l987) and the P-a correla- 
tion of lBlundell et al.1 (j 999) suggests sources with higher jet power 
have higher p. Also, 7 m in, often assumed to be 1 in previous mod- 
els, may in fact be higher and influence the luminosities, an issue 
we explore in § [Vj The minimum injected Lorentz factor 7 m in has 
taken various values f rom 1 to 10 4 for describing individual sources 
( lBlundelletal]|2006l) . 
We also assume 



u B (t) = 



rpi{t) 



(IWXr+l) 



and 



u e (t) = us{t)/r, 



(19) 



(20) 



that is, the ratio of the energy density in the particles to that in the 
magnetic field is a constant r = (1 + p)/4 , based on minimum 
energy arguments, adopted from lKaiser et alj dl997l) . 
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We now describe the determination of 71 and d^i/dj. First 



Vi dt 



t-tj+C/Vk 



t<tj 
t >tj 



(21) 



where 01 = 1 - (4 + /8)/[T (5 - /3)] and a 2 = 2/(3^ + 2-/3). 
Equation (O may be integrated to yield: 



A(t,ti) 



where 



4 (Tt 



A(t,ti) = i— [ (u B (t') + u c )F{t')dt' 
3 m e c L 



and 



F(t) = 



Fi(i) = t" ai/3 



*<*j 



F 2 (t) = (t - tj + Cj Vk) ~ a2 t>t s 



(22) 



(23) 



(24) 



The integral in A{t,t%) (equation \Ti\ ) has an analytic form as the 
integrand is proportional to a power of t for t < tj and a power of 
(t - tj + C/vT) for t > tj. Thus for ti < t ^ tj 



7i = 



and for ti ^ tj ^ t, 



with 



Fi(t)- 7 i4(t,t0 
Fl(tj)-7tj^(*j,*0 



(25) 



(26) 



(27) 



F a (t)-7A(t,*i) 

Then dji/dj can be found by differentiating equations l !25t and 
(HI. 

The total radio power at a frequency v can be then found by 
an integration over time: 



/•min[t,tj] 

P„(t) = / dP„ 
Jo 



(28) 



where dP„ is given by equation Q}. 



2.4 Inverse-Compton X-ray luminosity from the lobe 

We calculate the IC emission in the X-ray as a function of time 
as in An electron with a Lorentz factor of 7 boosts a 

CMB photon with frequency ccmb to frequency 7 2 ^cmb- We can 
make a simplification that all CMB photons are at the peak CMB 
frequency 5.879 x 10 10 Hz K" 1 x 2.7 K (1 + z). Then, the IC 
power at a frequency v is estimated as 

min[t,tj] I 3 

— a T oic- nMdV. (29) 
6n v 



p„(t) 



This is a simplification a nd we are n ot considering the full spec- 
trum of CMB photons but Nath ( 2010) shows that it is a reasonable 
measure of the X-ray luminosity, correct to within an accuracy of 
60 per cent. 

2.5 Evolutionary tracks 

We plot the time evolution of the length of the lobe in Figure [T] 
After time tj the length expands slower than when the jet was on. 
Power-linear size (P-D) diagrams for various source parameters 



a 
C5 



60 
O 
— 




Figure 1. Growth in the length of a lobe as a function of time. Qj 



1.5 X 



10 40 W, tj = 10 8 yr (top/red), Qj = 1.5 X 10 38 W, tj = 10 s yr 
(middle/green), Qj = 1.5 X 10 38 W, tj = 5 X 10 7 yr (bottom/purple). 
The growth rate of the lobe length changes once the jets turn off. 



are presented in Figure [2] The figur e is us e ful for comparing our 
mo del to the mo dels of [Kaiser et al] d 19971) . iBlundell et ail d 19991) 
and lNathl ( l2010l) which all present P-D diagrams. For the goals of 
this paper we simply need, however, the time evolution of the X- 
ray and radio luminosities. We will consider the X-ray luminosity 
at 1 keV and the radio luminosity at 151 MHz. The time evolution 
at these two frequencies is presented in Figure[3] 

In Figure[5]we see that while the jet is switched on, the radio 
luminosity declines and the X-ray luminosity rises. After the time tj 
when the jet turns off the radio luminosity precipitates very quickly 
and the X-ray luminosity falls but at a less rapid rate than this. At 
higher redshifts, the X-ray luminosity is stronger due to the (1 + z) 4 
rise in the energy density of the CMB. The luminosities also fall 
more rapidly at higher redshifts due to increased IC losses from 
the CMB. In Figure|4]we present the length of time a radio source 
continues to radiate in the X-ray after jet activity has ceased. At 
z = 2, a source will radiate in the X-ray for about 30 Myr after the 
jet activity stops at tj — 10 s yr, since there is no new injection of 
particles and existing particles had their Lorentz factors reduced to 
less than 7 ~ 10 3 . 



2.6 Model parameters 

Model parameter case [A] has tj = 10 s yr, p = 2.14, 7 m in = 1, 
po = 1.67 X 10 -23 kg m~ 3 and j3 = 1.5. These parameters are in- 
ferred from observations and similar to parameters used in previous 
models for radio sources. Parameters such as the minimum injected 
Lorentz factor, 7 m i n = 1, however, are no t well constrained b y 
observations and may range from 1 to 10 dBlundell et ai] |2006). 
Other sets of model parameters we investigate (§ |7J are: [B] tj — s> 
5 x 10 7 yr, [C] tj -> 5 x 10 8 yr, [D] p -»• 3, [E] 7min -> 2000, 
[F] p — )• 3 and 7 min -J- 2000, [G] p correlated with Qj (see § 0, 
[H] po -> 1.67 x 10 -24 kg m -3 , [I] /3 -> 2 [J] birth function (see 



© 2010 RAS, MNRAS OOO.rjJfJJ] 



Inverse-Compton ghosts and giant radio sources 5 




2.0 2.5 

Log 10 D (kpc) 

Figure 2. Power-linear size (P-D) diagram for the evolution of the source 
in the X-ray (1 keV) and the radio (151 MHz). Radio luminosities are 
dashed, X-ray luminosities are solid. Tick marks represent time intervals 
of 0.2ij increments, starting from 0.2ij. [A/blue] Qj 



■ 10 s yr, z = 
1, [C/red] Qj 



, [B/green] Qj = 
1.5 x 10 40 W, 



1.5 x 10 38 W, 
1.5 x 10 38 W, fc = 10 s yr, 



Qj = 1.5 x 10 38 W, tj = 5 x 10 7 yr, z - 



10 s yr, 



2, [D/purple] 



so 
c 




Figure 4. The length of time a source continues to radiate in the X-ray 
(1 keV) after the jet has turned off as a function of redshift (similar for all 
jet energies). Jet lifetimes of tj = 5 X 10 7 yr, 10 8 yr and 5 x 10 8 yr 
(bottom, middle, top / green, blue, purple) are used for the figure. 



-3.0, and [L] 



§ l3.2t Az — > 1, [K] jet-power power-law —2.6 
jet-power power-law —2.6 — > —2.0. 



5/J 
O 




100 150 

t (Myr) 

Figure 3. Power-time (P-t) diagram for the evolution of the source in the 
X-ray (1 keV) and the radio (151 MHz). Radio luminosities are dashed, 
X-ray luminosities are solid. [A/blue] Qj = 1.5 X 10 38 W, tj = 10 s yr, 
2 = 2, [B/green] Qj = 
Qj = 1.5 xlO 40 W,tj 
tj = 5 X 10 7 yr, z = 2 



1.5 X 10 38 W, tj = 10 8 yr, 2 = 1, [C/red] 
= 10 8 yr, 2 = 2, [D/purple] Qj = 1.5 X 10 38 W, 



2.7 Comparison with other models in literature 

Relevant earlier work in the subject goes back to 
Gopal - Krish na" et alj (j989), which produced the first quanti- 
tative study of the evolution of old radio lobes and their fading via 
IC processes. A number of models since have been developed for 

the evolution of the size and lum inosity of radio lobes . 

The m odel p resented in iKaiser & Alexander! d 19971) and 
iKaiser et al.1 d 19971) provides an analytical description of the evo- 
lution of FR II radio sources that includes adiabatic, synchrotron, 
and IC losses on the source's luminosity. The model assumes a self- 
similar expansion of the jet and lobes, with a constant axial ratio 
throughout evolution. The evolution does not go beyond the life- 
time of the jet. 

1999) 



To investigate spectral dependencies Blundell et al 



created a model ba sed on that of IKaiser & Alexander 1 1997 ) and 
IKaiser et all dl997l) . however, their model considers a more com- 
plex injection into the lobe governed by breaks in the energy dis- 
tribution determined by the dwell time and magnetic field strengths 
in the hotspot, rather than assuming a constant injection index. The 
model also considers the hotspot pressure rather than the lobe pres- 
sure to govern adiabatic losses out of the hotspot in order to find 
the lobe luminosity. Axial ratios in this model increase with time 
(which correctly accounts for the observed larger increase for more 
powerful jets), as the hotspot pressure does not scale in the same 
way as the lobe pressu re. This model bet ter r eproduces observa- 
tions o f jets, such as in iLeahv et alj dl989h and lLeahv & Williams! 
( 1984), which found that axial ratios are larger in sources with 

high er jet powers. 

iManolakou & Kirld d2002l) extend the model of lBlundell et al.l 
dl999F Instead of specifying the electron distribution that enters 
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the lobe, the model accounts for electrons accelerating by the first- 
order Fermi process at the termination shock and then propagating 
through the hotspot into the lobe. To match P-D observations, a 
reacceleration process during propagation through the head is in- 
cluded which compensates for the adiabatic losses between the ter- 

min ation shock and the lo be. 

iBarai & Wiital d2006r) andlWang & Raised d2008h compare the 
thr ee models oflKaiser & Alexander! jl997l) . lBlundell et al.ld~1999l) . 



an d lBlundell etaLH 1999). Monte Carlo simulations for each model 
are used to predict radio powers, sizes, redshifts and spectral in- 
dices of an artificial sample, which are compared to the data of the 
low-fr equency radio survey 3CRR, 6CE and 7CRS. Barai & Wiita 
find that no existing model can give acceptable fits to 
all the properties of the surveys considered, and the simplest 
iKaiser & Alexander! J 1997b model is somewhat better at fitting the 
data. IWang & Kaise evolve some of the properties of FR 
II sources with redshift t o have the artificial samples to fit the ob- 
servations, and also finds IKaiser & Alexander! d 19970 best matches 
observations. 

The model of lNathl feoich is the first to go beyond the stop- 
ping of the jet after a time t$ and consider the X-ray emission of 
the lobes from IC scattering of the C MB photons. The mod el is a 
variant of lKaiser & Alexander! d 19971) and lKaiser et alj d 19970 . also 
assuming a constant axial ratio. The model does not evolve simply 
according to the self-similar evolution of lobe length determined by 
the characteristic length scale. 

The model in the prese nt work also goes beyond the cessation 
of the jet in the formalism of Nath (2010) and considers the growth 
of the lobe as described bv lBlundell eTal] jl999h where a compact 
hotspot with a pressure distinct from the lobe pressure that deter- 
mines lobe length growth feeds the lobes, and in which axial ratios 
grow over time during jet activity. The model is also analytic rather 
than numeric. We are able to go beyond the self-similarity expan- 
sion after the jet stops (i.e. no longer rely on a characteristic length 
scale). We are not required to fix a constant axial ratio during jet ac- 
tiv ity and solve the pairs of differential equations as in equation (6) 
of iNathl d201dh to solve for the evolution of length and pressure, 
which results in the volume evaluated as Vi(t) = TvLj (t) 3 / (AR 2 ) 
not precisely agreeing with the volume of the l obe determi ned by 
integration of volume elements (equation (8) of lNathl d20100 ). 



3 OBSERVATIONALLY INFERRED DISTRIBUTION 
PARAMETERS 

In this section we discuss the radio luminosity function (RLF) of 
Fanaroff and Riley Class II radio sources, the empirically-inferred 
birth function of such sources and the empirically-inferred distri- 
bution of jet energies which are necessary to find the distribution 
of the IC ghosts and the completely dark, i.e. non-observable, ra- 
dio and X-ray lobes. An IC ghost refers to a galaxy that has its jets 
turned off but still emits X-ray radiation due to the upscattering of 
CMB photons. An observable IC ghost is the term we use to call 
an IC ghost that radiates above a given X-ray flux limit (but lacks 
detectable synchrotron radiation). In § 13161 we use the normative 
model parameters of case [A]. In §[7] we explore model parameter 
cases [B]-[L]. 

3.1 RLF of FR II sources 

We will u se the RLF of high radio power FR II sources deter- 
mined by IWillott et al.1 feOOlh from the 7CRS, 6CE and 3CRR 



samples, to aid in predicting the density of radio lobes including 
IC ghosts and completely non-observable ones. The radio luminos- 
ity function describes the space density per unit comoving volume 
of sources as a function of luminosity as derived from the sur- 
veys. The luminosity function was converted from the cosmology 
H = 50 km s" 1 Mpc~\ fi M = 0, f2 A = , » k = 1 to t he 
cosmology we use in this paper with the relation jPeacocklll985h 



In \ dVl In \ dV ' z 



(30) 



where the indices refer to a cosmological model, Pi is the luminos- 
ity derived from the flux density and redshift in model i and Vi is 

comoving volume; 

The RLF of IWillott etakl feOOll) will reflect only those radio 
galaxies that are abov e the flux limit (at low radio frequencies ~ 
151 MHz) of 0.5 Jy dGopal-Krishna & Wiitall200lh . which is the 
lowest flux limit of the surveys they used to determine the RLF. 



3.2 Birth function of radio galaxies 

We assu me the empirically-i nferred birth function of radio sources 
given bv lBlundell etai1dl9990 : 



p(z)dz oc e 



dz 



(31) 



with zi = 2.2 and Az — 0.6. Most sources are born during the 
quasar era, z — 1.5-3. By a "birth" w e mean jet activity in itiates. 
The birth function was determined by iBlundell et alj d 19991) such 
that simulations of sources best match 3C and 7C data. 

It is useful to convert the birth function to the probability den- 
sity function (PDF) for the age of radio sources as a function of 
redshift. This is accomplished through the equation 



PA(t,z)dt - 



p(z t ) 



p{z') dz' 



-dz 



(32) 



where pA(t, z) gives the probability density that a source is age t 
at redshift z and zt is the redshift of a source at its birth for it to 
appear age t when observed at a redshift of z. The PDF for the age 
of radio sources at z — 0, 1, 2, 3, 4 is presented in Figure|5] 



3.3 Distribution of jet powers 

We use the empirically-inferred p robability density function of jet 



energies bv lBlundell et 



-interred pi 
al. ( 1999). 



The distribution is given by 



p(Q)dQ oc Q- 



(33) 



. The lower- 



withQ min = 5xl0 37 W < Q < 5xl0 42 W = Q 
energy jets are most abundant according to this power law. At high 
redshifts we may be only seeing the brighter sources (those that are 
younger and/or have higher jet power) above the radio flux limit, 
meaning that a large number of radio sources may go undetected. 
We will also c onsider modifying the power-law index from —2.6 
to —2 and — 3. IWang & Kaiserl d2008) find that a power-law index 
of —2 better matche s observed P — D — z and a distributions. 
Bara i & Wiital (2006) on the other hand find that a steeper power- 
law index of —3 is preferable. 
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Figure 5. Probability density function (PDF) for the age of radio sources at 
various redshifts. z = 0, 1, 2, 3, 4 are shown. The PDF of ages is useful for 
determining the density of IC ghosts as well as completely dark X-ray and 
radio lobes. 



4 OBSERVABLE DISTRIBUTION OF 

INVERSE-COMPTON GHOSTS AND ACTIVE 
SOURCES 

We are interested in the distribution of all radio galaxies that have 
ongoing jet activity, the IC ghosts and galaxies with dead radio 
lobes that no longer exhibit detectable radiation. Dead lobes refer 
to lobes that are completely X-ray and radio dark. Although unde- 
tectable, the fossil radio lo bes are still an importan t comp onent of 
the intergalactic medium: lEnfilin & Gopal-Krishnal fcOOlh demon- 
strated that radio plasma in the lobes of galaxies, after jet activity 
has ceased, can be revived even up to 2 Gyr later by compression 
in a shock wave produced by large-scale structure formation. The 
major steps of the analysis are summarized in a flowchart in Fig- 
ure [6] 

The ratio of IC ghosts to radio sources with jets on at a given 
redshift z is: 



IC ghosts 



p\(t, z) dt 



sources with jets on 



p A (t,z)dt 



(34) 



where tj is the jet lifetime and t max is the time at which the X-ray 
luminosity goes to 0, or, if we care about the observable distribution 
of the ghosts, then the time at which the X-ray luminosity drops 
below some flux limit. We may alter the integration limits in the 
numerator accordingly to derive the ratio of sources with dead lobes 
to sources with jet activity (specifically, integrate from t max to oo). 

Thus, if we have an accurate idea of the total density of radio 
sources with current jet activity at each redshift, then we can esti- 
mate the density of IC ghosts and dead lobes as well using a ratio 
as in equation d34b . with appropriate integration limits. 

We can use the RLF of FR II objects described in § 13.11 to 
find the total density of radio sources with jets on, p s j (z). Recall 
that the RLF will only account for those radio galaxies that are 



Use empirically inferred birth function and distribution 
of jet energies consistent with the RLF of FR II sources 



From the birth function find the prob- 
ability density function (PDF) of 
sources with different ages given redshift 



Use model developed in the paper to predict the age 
at which a source falls below radio flux limit used to 
determine the RLF, given a redshift and jet power 



Infer the total population of FR II sources from the 
PDF of ages and the fraction accounted for by the RLF 



At this point we can consider any subpop- 
ulation (for example active sources, ghosts) 
by specifying objects within a range of ages 



Calculate: (1) number of sources observable per 

square degree given a flux limit, (2) contribu- 
tion to diffuse X-ray background, (3) volume fill- 
ing factor of lobes, (4) X-ray luminosity function 



Figure 6. A flowchart summarizing the major steps in the analysis of § [5] 
and|4] 



above the flux limit of the radio survey. In other words, because the 
radio luminosity of a source decreases with time in our model, the 
RLF only accounts for galaxies that are younger than some time 
t'(Qj, z) ^ tj at each redshift in the evolution of the radio source 
with jet power Qj and redshift z. Then, to find the total density of 
radio galaxies with jet activity at a given redshift z, the ratio 



sources with jets on 
FR II from RLF 



PA(t)dt 



/ Q " maX p(Q) / 

jQj min JO 



t'(Q,z) 



p A (t)dtdQ 



(35) 

is used. Here we are averaging over the empirically-inferred jet 
powers. The density of FR II objects accounted for by the RLF 
is found by a simple integration of the RLF over luminosity. Mul- 
tiplying this density by the ratio in equation d 3 5 1 > gives p s j (z) that 
was sought. Then, at any redshift, the density of all objects active 
or dead, p to t(z), is: 



ptot(z) = 



p A (t,z)dt 



(36) 



An X-ray limiting flux of 3.0 x 1CT 1J J iri s _1 will be 
considered, similar to that of Chandra X-ray Deep Field North 
(CDFN) survey. Also, a radio limiting flux density of 7.55 x 



10" 



J m s will be considered, corresponding to the flux 



limit of the 7CRS survey which had the lowest flux limit of the 
surveys used in the determination of the RLF of FR II sources. 
We are interested in the number of X-ray ghosts above the X- 
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ray flux limit per square degree. This is obtained by the equation: 



N = Utt 



Pghosts>x.f .1. 

dV c 



(37) 



where p g hosts>x.f.i. is the comoving number density of the IC 
ghosts above the X-ray flux limit obtained from multiplying p s j 
and the ratio in equation i34i . 



dV c = 4tt- 



£l(z) 2 



-dz 



(38) 



Ho (1 + z) 2 £(z) 

is the comoving volume of a shell at z of thickness dz, where Dl 
is luminosity distance and 

E(z) = (n M (l + zf + Q.a + n k (l + z) 2 ) 1/2 . (39) 

(see lCappil fefJOlh ). We integrate from z = to z = 4 to find the 
number density. 

We predict ~ 29 powerful double-lobed galaxies per square 
degree visible in the X-ray above the X-ray flux limit (of the 
CDFN survey), of which ~ 4 (approximately 13 per cent) are IC 
ghosts. Our prediction is 17 per cent of the number d ensity of ex- 
tended X-ray sources estimated by the CDFN survey: iBauer et all 
(2002) find a surface density of IQT^Jj (1 a) extended X-ray 
sources (including X-ray clusters) at a limiting soft-band flux of 
3 x 10~ 16 erg cm" 2 s _1 = 3.0 x 10 -19 J m~ 2 s _1 inferred from 
the ~ 1 Ms CDFN observation. 

We also calculate the expected number density of sources 
for an X-ray flux limit of 2 x 10 -15 erg cm -2 s _1 = 2 x 
10~ 18 J m~ 2 s -1 , which was the depth of the Sub aru-XMM Deep 
Field North observation. Finog uenov et alj ^20 id) identify 57 X- 
ray cluster candidates and 6 sources possibly due to X-ray emission 
of radio lobes in the Subaru-XMM Deep Field which surveyed a 
1.3 square degree region (about 5 sources per square degree). Our 
model predicts 1 X-ray source per square degree above the flux 
limit. The predicted number of observable IC ghosts and total radio 
lobes per square degree in the X-ray as a function of the X-ray flux 
limit is presented in Figure[7] 

Also of interest are the luminosity ratios of t he X-ray (1 keV) 
to rad io (151 MHz) predicted by our model. ICelotti &" Fabian 
( 2004) use this ratio to convert the RLF of FR II objects to an X-ray 
lu minosi t y func tion (the ratio is taken to be a conservative value of 
1). iNathl ( 1201 01) develops a model for the evolution of the X-ray 
and radio luminosities of double-lobed radio galaxies and takes the 
time-averaged ratio over the duration in which the sources were 
both X-ray and radio bright (in the rest frame) to also obtain an 
XLF for extended X-ray sources. We plot the ratios at t — 0.05tj 
and t = ij of the jet as a function of redshift in Figure[8]for various 
jet energies. Fo r comparison, the time-avera ged ratio 0. 14(l + z) 38 
of INathl j2O10b is included as well, which iNath] d2010l) obtains by 
averaging the ratio of luminosities in X-ray and radio frequencies 
of that paper's model over a time until the X-ray or radio emis- 
sion drops rapidly. We also plot ra tios obtained from real obser- 
vations from 4C 6 0.07 (z = 3.79) dSmail et al.ll2009h . PKS 1138- 
262 (z = 2.156) dCarilli et alj | 2002h . PKS 0156-252 (z = 2.09), 
PKS 0406-244 (z = 2.44), PKS 2036-25 4 (z = 2.00), PKS 2048- 
272 (z = 2.06) dOverzier et ail 120051) . 3C 432 (z = 1-785) , 
3C294 (z = 1.779), 3C 19 1 (z = 1.956) jBriund et alj|2006t) . 
6C 0905+39 (z = 1.883) ferlund et all 120081), MRC2216-20 6 
(z = 1.148), MRC 09 47-249 (z = 0-854 ) dLaskar et alj|20ich . 
4C 23.56 (z = 2 4 8) d Johnson et ai]|2007l) . 4 C41.17 (z = 3.8 ) 
dScharf etalj|2003l) and 3C98 (z = 0.0306) dlsobe et alj|20051) . 
The extended X-ray lobe luminosities (1 keV rest frame) are ob- 
tained from the mentioned papers (recalculated to our cosmology) 




-20 -19 -18 

Log |f) X-ray flux limit (J m~ 2 s -1 ) 



-17 



Figure 7. The number per square degree of all radio lobes visible in the 
X-ray (top/blue) and the number of observable IC ghosts (bottom/green) as 
a function of X-ray flux limit. 



and the radio luminosities are estimated from NASA/IPAC Extra- 
galactic Database (NED) photometric data. The size of the sample 
is 15, but we can already see large scatter in the ratios. Our model 
is also able to predict a wide range of ratios. Some sources how- 
ever appear to be above the ranges predicted by our model. We are 
only showing the range of ratios while the jet is on. Once the jet 
turns off the radio luminosity plummets very rapidly and the X-ray 
to radio ratio tends towards infinity. It may be possible that some 
of the sources with high ratios have their jets turned off recently. 
Or, it may be possible that the sources with high X-ray to radio 
luminosity ratios have jet lifetimes greater than 10 s yr or that envi- 
ronment or injection spectrum parameters of these sources deviate 
significantly from those of case [A]. We can obtain higher ratios if 
we increase p or lower th e surrounding density. Also, if /3 — 2, the 
value used in INathl d2010l) . the ratios are found to not depend on jet 
power. We are likely underestimating the X-ray to radio luminosity 
ratios of the 15 sources because the NED radio luminosity includes 
the lobes, hotspot and nucleus rather than just the lobes. 

4.1 Diffuse X-ray background contribution 

The diffuse X-ray background contribution from accumulated 
Compton lobes and ghosts is also calculated. It is important to 
check that this flux does not exceed the observed CXB. The diffuse 
X-ray flux (in J m~ 2 s _1 deg -2 ) can be found by an integration 
across comoving volume shells and hence z. 



F = 



f(t,z)p tot {z)dV c (40) 



with 



/(*,*) 



2L x (t) 
4ttA 2 (z) 



x p A (t, z) dt, 



(41) 



for example, gives the X-ray flux of the IC ghosts. The luminos- 
ity will depend on jet power Qj, so we bin jet power, calculate F 
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Figure 8. X-ray (1 keV) to radio (151 MHz) luminosity ratios as a function 
of redshift. A jet lifetime of tj = 10 s yr is used for the figure. The ratios in 
our model at t = 0.05fj (bottom set of lines) and t = t j (top set) are shown 
for [1] Qj = 1.5 X 10 38 W (dotted/red), [2] Qj = 1.5 X 10 40 W (dot- 
dashed/green), [3] Qj = 1.5 X 10 42 W (dashed/purple). Once the jet turns 
off, radio luminosity drops rapidly and the ratio tends to infinity. The time- 
averae ed ratio (0.14(l + z) 3 - 8 ) over the duration in which the sources were 
both X-ray and radio bright in the rest frame calculated in Nath (2010) is 
also shown (solid/blue). The blue points represent actual sources (see §0. 

at each and take a weighted average according to the empirically- 
inferred distribution of jet powers. We calculate luminosities at 
1 keV observed energies. 

From studies of the CXB bv lHickox & Markevitchl d2007t) and 
Iffickox & Markevitchl d2006t) . the flux density of the CXB (after all 
resolved X-ray point and extended sources have been removed) at 
~ 1 keV is approximately 4 X 10 -12 erg cm -2 s -1 deg -2 = 
4 X 10 -15 J m -2 s -1 deg -2 . The X-ray flux from all the sources 
predicted by our model is 9.7 x 10 -17 J m -2 s -1 deg -2 , about 
2.5 per cent of the observed but unresolved CXB. Of course, many 
of the objects in the calculation of X-ray flux are above the flux 
limit and should be excluded. The flux from all extended sources 
below the X-ray flux limit is 7.3 x 10 -17 J m -2 s -1 deg -2 . Comp- 
ton ghosts contribute about a tenth of the total X-ray flux: 1.1 x 
10 -17 J m -2 s -1 deg -2 , of which 3.9 x 10 -18 J m -2 s -1 deg -2 
comes from sources that are below the X-ray flux limit. 

The diffuse X-ray background contribution from lobes of 
powerful double radio sources is found to make up 0.02 of the unre- 
solved CXB. IC scattering of the CMB off the nucleus and hotspots 
of sources as well as X-ray clusters and FR I objects also contribute 
to the CXB. 



5 VOLUME FILLING FACTOR OF THE UNIVERSE BY 
LOBES 

Next the volume filling factor £ of all lobes and relic lobes is cal- 
culated at each redshift. 

C= / 2Vi(t)p A (t,z)dtX (1 + z) 3 Ptot (z) (42) 
Jo 



gives the volume filling factor at z for sources with age less than 
T (again, the limits of integration altered appropriately determine 
which objects we care about). T is set to the age of the universe at 
redshift z if we care about all objects, and is altered appropriately 
to account for just active lobes or IC ghosts. We allow for expan- 
sion of the volume V\(t) only until the time the lobe pressure comes 
into equilibrium with the ambient pressure, calculated from the am- 
bient density profile assuming a temperat ure of 10 5 K (temper a- 
ture of the warm-hot intergalactic medium dCen & Ostriker 2006)), 
at which point the volume remains constant. The ratio of the lobe 
pressure at time t to the ambient pressure at the distance of the end 
of the lobe at time t stays well above 1 for times beyond which 
the lobe has Mpc scale length. The lobe volume and source den- 
sity will depend on jet power Qj, so in order to find £ we bin jet 
power, calculate £ at each and take a weighted average according 
to the empirically-inferred distribution of jet powers. The factor of 
(1 + z) 3 in the expression for £ comes from converting comoving 
volume to proper volume. 

We also consider the relevant volume the lobes are situated in, 
namely the volume of filamentary structures of the universe. We 
divide the volume filling factor by lobes by the vo lume fraction of 
the wa rm-hot intergalactic medium in figure 1 of ICen & Ostriked 
(2006) to account for the relevant volume. Figure [9] shows the vol- 
ume filling factors as a function of z. At z = 2 and z — 3, we 
predict volume filling factors of all lobes around 0.03 and 0.02 re- 
spectively. For a longer jet lifetime, tj = 5 x 10 8 yr, the predicted 
volume filling factors at these redshifts are 0.3 and 0.2. While lobes 
may not significantly fill the entire universe, it is possible that they 
fill a large fraction of the filaments in which they are located. 

iGopal-Krishna & Wiital feOOll) give a simple estimate for the 
volume filling factor by lobes as well. Their paper estimates £ = 
0.01 for tj = 10 8 yr and C = 0.53 for tj = 5 x 10 8 yr (at 
z = 2). IGopal-Krishna & Wiital l l200ll) do not use a birth func- 
tion of radio sources and estimate the number of dead lobes from 
tq e /tj, where t qe is the duration of the quasar era. This estimate 
gives slightly higher ratios at z = 2 than estimated by the birth 
function due to the simpl er age distribution of galaxies . Typical 
lobe volumes are larger in IGopal-Krishna & Wiit j J200lh than es- 
timated by our model since their paper assumes a constant ax- 
ial ratio of 2.5 (according to this paper's definition of axial ra- 
tio), i.e. assuming more spherical sources. The question of vol- 
ume filling factor of lobes is al so ad dressed with more care in the 
works by |Barai& Wiital d2007l) and lBarail d2008h . with which the 
results of our work agree. The former paper investigates the vol- 
ume fi ll ing factor by lobes usi ng th e models o f Kaiser & Ale xander 
dl997l) . iBlundell et alj Jl999l) and iManolakou & Kirk! d2002l) . and 
also modifications by incorporating a variable hot spot size grow- 
ing with the source age. Table 7 of their paper gives the relevant 
volume fraction results for the models with a number of parame- 
ters varied. A wide range of relevant volume filling factors can be 
produced by modification of the parameters, with the cumulative 
relevant volume filling factor of radio galaxies over the quasar era 
around 0.05. The later paper incorporates radio lobe growth into a 
numerical cosmological evolution, and finds a volume filling factor 
of 0.10-0.30. 

The total number of particles placed into the IGM because of 
the formation of the lobes is 

AWtidcs = J o ' (dV(ti) J*™" n (ti) 7 ~ P *y) («) 

N is found to be between 10 65 to 10 67 particles for the different jet 
energies. The density of these particles in the lobes at the time the 
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Figure 9. Volume filling factors of sources as a function of redshift. The 
solid lines from top to bottom refer to: [1/green] all radio galaxies including 
completely X-ray and radio dark ones, [2/blue] IC ghosts and sources with 
radio jets on, [3/red] IC ghosts. The corresponding dashed lines show the 
volume filling factors if we consider the relevant volume of the filamentary 
structures of the universe. 



jet turns off is on the order of 10~ 32 -10~ 31 kg m" 3 . Therefore, the 
ratio of the density of particles in the lobes to the critical density 
of the universe f2 ra dio lobes is at least 6 orders of magnitude below 
unity, even for a lobe volume filling factor of 1. 



6 X-RAY LUMINOSITY FUNCTION OF EXTENDED 
X-RAY SOURCES 

The X-ray luminosity function for extended X-ray sources may be 
directly found from our model since we have described the density 
and the age and jet energy distributions of the sources and have a 
model for the evolution of the lu minos i ty. W e are not required to 
take a time-average of ratios as in iNathl ( 120101) to convert the RLF 
to an XLF. We present in Figure [10] the XLF for all extended X- 
ray sources derived from our model and syn thetic population. F or 
comparison, the XLF for X-ray clusters from lMullis et al.l {2004) is 
also plotted. The evolution of the cluster XLF is not well known at 
luminosities < few x 10 44 erg s _1 and beyond redshift z — 0.8, 
however, observed volume densities for 0.6 < z < 0.8 and for 
luminosities above few x 10 44 erg s -1 are signifi cantly lower than 
those at the present population (Mul lis et alj2004l) . 

The XLF derived from our model, when integrated, gives 
a higher number density than does integrating the original RLF 
because it accounts for the predicted total number of luminous 
sources, whereas, we have mentioned earlier, the RLF has no 
knowledge of the evolution of radio sources and therefore when 
integrated only accounts for the density of the sources above the 
lowest flux limit of the samples used in its construction. 
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Figure 10. The radio (151 MHz) and X-ray (1 keV) luminosity function 
for various sources. The bottom, middle, top dashed (blue) lines refer to the 
RLF of FR II objects derived from the 7CRS, 6CE and 3CRR surveys at 
z = 0.2, 1, 2 respectively. The bottom, middle, top thick solid (red) lines 
refer to the XLF of all extended X-ray sources derived from our model 
(tj = 10 8 yr) and age and density distribution at z = 0.2, 1, 2 respec- 
tively. Dotted thin lines show the IC ghost part of the population and solid 
thin lines shows the cu rrently active galax y population. For comparison, the 
XLF of X-ray clusters I Mulli s etall2004l) (dot-dashed/grey) is shown (top, 
middle, bottom correspond to z = 0.2. 1, 2 respectively). 



7 VARYING MODEL PARAMETERS 

We test how sensitive our predictions are to changes in the model 
parameters. Varying the parameters can indeed have a significant 
effect on the distribution of double-lobed sources and IC ghosts, 
since different model parameters change the times sources go be- 
low the radio flux limit, and density parameters and jet lifetimes 
also increase or decrease the volumes of the sources. The parame- 
ters we alter are the jet lifetime tj, the injection index p, the mini- 
mum energy of injected electrons 7 m in and parameters describing 
the surrounding AGN environment. In addition, we also investi- 
gate the effect of the birth function. We observe how our results 
change with shorter and longer jet lifetimes of tj = 5 x 10 7 yr 
and tj — 5 x 10 s yr (cases [B] and [C]). The injection index 
p — 2.14 we considered is not significantly steep so we also 
consider the steep case ([ D]) of p — 3 (observational work by 
lAlexander & Leahvl < fl987h find 2 sC p sC 3). Additionally, we con- 
sider correlating p with jet power Q\ (c ase [G]), inspired by the 
P—a correlation of lBlundell et alj dl999l) . A simple linear correla- 
tion between p and log 10 Qj is assumed, withp = 2 corresponding 
to Q m in and p — 3 to Q max - We consider injected electrons to 
have 7 m in = 2000 (cases [E] and [F]), greater than the 7 ~ 10 3 
factors required for the X-ray radiation, in addition to the original 
case ([A]) which had j m i n = 1. We also examine the importance 
of the surrounding density and the effect of reducing the density 
by an order of magnitude by setting po = 1.67 x 10~ 24 kg m~ 3 
(case [H]), as well as considering an ambient density profile given 
by a steeper power law with f} — 2 (case [I]). We consider a wider 
birth function with Az = 1 (case [J]). Finally, we consider modi- 
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Figure 11. Volume filling factors of the filamentary structure of the universe 
by lobes at redshifts z = 1, 2, 3 for various model parameters. Case [A] 
original parameters investigated (tj = 10 s yr, p = 2.14, 7 m i n = 1, 
po = 1.67 x 1(T 23 kg m~ 3 , = 1.5), [B] t j -> 5 x 10 7 yr, [C] tj -> 
5 X 10 8 yr, [D] p 3, [E] 7 min -> 2000, [F] p -5- 3 and 7 min -> 2000, 
[G] p correlated with Qj, [H] p -» 1-67 X 10~ 24 kg m" 3 , [I] /3 -> 2, 
[J] birth function Az — > 1, [K] jet-power power-law —2.6 — > —3.0, [L] 
jet-power power-law —2.6 — > —2.0. The horizontal dashed line shows a 
volume filling factor of unity. In some instances the predicted volume filling 
factor exceeded unity, which means that the ejected material from different 
sources overlap. 



fying the power-law distribution of jet energies to —3.0 and —2.0 
(cases [K] and [L]). In Figures [TT1 II 21 [T3l we present how the vol- 
ume filling factor, source number density above the X-ray flux limit 
and the unresolved X-ray emission change with varying the model 
parameters. 

In the original case [A] we predict ~ 29 lobed-galaxies per 
square degree visible in the X-ray, of which ~ 4 are IC ghosts, 
above the X-ray flux limit of the CDFN survey. The radio sources 
contribute about 2 per cent of unresolved CXB and the volume fill- 
ing factor of the lobes is ~ 0.03-0.02 at z = 2-3. 

The jet lifetime is an important parameter that affects the 
predicted distribution of the sources. A larger jet lifetime means 
that sources will grow longer while there are also fewer predicted 
sources with jets turned off during the quasar era. The predicted 
volume filling factor exceeds 0.25 at redshifts of the quasar era 
for an increased jet lifetime of 5 x 10 8 yr, while the volume fill- 
ing factor decreases to £ ^ 0.01 for a decreased jet life time of 
5 x 10 7 yr. Evidence that a significant fraction of sources have 
long jet lifetimes on the order of 5 x 10 s yr will mean that double- 
lobed sources can account for 20 per cent of the unresolved CXB 
and 60 per cent of the extended sources in the CDFN survey. 

Increasing the injection index to p = 3 has the effect of re- 
ducing the radio luminosities and consequently sources are more 
likely to be below the radio flux limit and hence a larger number 
of total sources in the underlying distribution are predicted. Con- 
sequently, the predicted volume filling factor of sources at redshift 
2 are significantly higher, surpassing unity in the quasar era and 
even at z — 1. The contribution of lobes to the unresolved CXB 
is 60 per cent (Figure [T3t. Despite the large volume fraction, the 
sources are less bright in the X-ray compared with [A] and the num- 
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Figure 12. Number of sources predicted to be visible above X-ray flux limit 
of CDFN survey for various model parameters, 'a' stands for all sources 
(switched on and switched off), 'g' stands for IC ghosts. [A]-[L] same as in 
Figure fTTI The horizontal dashed line shows the number of extended X-ray 
sources expected to be observed above X-ray flux limit of CDFN survey. 
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Figure 13. Expected X-ray flux from radio sources expected for various 
model parameters. T stands for the total flux, 'u' stands for the unresolved 
part (i.e. due to objects below the X-ray flux limit). [A]-[L] same as in 
Figure rm The horizontal dashed line shows X-ray flux from the CXB after 
all resolved X-ray point and extended sources have been removed. 



ber of sources observed per square degree is 9 at the flux density 
limit of the CDFN survey. 

The minimum electron injection energy determined by 7 m i n , 
when increased to 7 = 2000, leads to higher X-ray and radio lumi- 
nosities than our original case, meaning there are fewer sources pre- 
dicted to go undetected and hence a smaller number of total sources 
in the underlying distribution and volume filling factor. The volume 
filling factor is still as high as 0.017 during the quasar era. Addi- 
tionally, setting 7 m in = 2000 results in the X-ray luminosity falling 
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Figure 14. The evolution of Lorentz factors since time of injection into the 
lobe. For a jet lifetime of tj = 10 8 yr, Qj = 10 38 W, z = 2 we show 
how a particular 7 evolves as a function of time. We plot the evolution of 
7min = 2000 when it is injected into the lobe at times O.ltj (purple), 0.5tj 
(green) and tj (blue). We also show the evolution of 7 max = 10 6 injected 
at time tj (blue, dashed), defining an upper limit of observable 7 factors 
after the jet stops. The vertical dashed line shows the time the jet switches 
off. 



less steeply for the majority of the time the source is an IC ghost, 
meaning that it is above the flux limit for a longer period of time 
and hence more observable. The ratio of observable IC ghosts to all 
observable sources for 7 m i n = 2000, at the flux limit of the CDFN 
survey is 30 per cent versus 13 per cent for 7 m i n = 1. 

By changing both p and 7 m i n appropriately it may be possi- 
ble to have a typically high injection index and predict reasonable 
values for the observable distribution (number density, CXB con- 
tribution, volume filling factor, especially after the quasar era) of 
double-lobed sources. Such is the case when p — 3 and Tmin = 
2000, where we are not underestimating radio luminosities as in 
the p = 3 and 7 m i n = 1 case. The minimum injected Lorentz fac- 
tor 7 m in, whose typical value is not well known, appears to be an 
important parameter for understanding the luminosity evolution of 
powerful double-lobed sources. Figure [T4l shows the evolution of 
Lorentz factors since the time of injection. 

In correlating injection index p with jet power Qj, we have 
made the least powerful sources brighter and the most powerful 
sources dimmer. This could, in principle, alter the volume filling 
factor and number density in either direction. For our correlation 
we described earlier in this section, the number density and vol- 
ume filling factor are increased somewhat from the original case, 
and double-lobed sources account for 27 per cent of the extended 
sources of the CDFN survey and 6 per cent of the unresolved CXB. 

The ambient density of galaxies and its evolution, which is 
still not well understood, could significantly affect the estimate for 
the volume filling factor of filaments by radio lobes. Decreasing the 
ambient density (reducing po by a factor of 10) allows for sources 
to grow larger and initially more luminous, but luminosity also falls 
more quickly due to increased adiabatic losses. Volume filling fac- 
tors are increased from the original case to above unity. The observ- 



able number density of sources is higher by a factor of 2, however, 
than observed by the CDFN survey. The unresolved contribution to 
the CXB is 45 per cent. We also investigate the effect of a steeper 
density profile with /3 = 2. The surrounding density is larger than 
the original case for distances r < ao — 10 kpc and smaller at 
distances r > ao. The lobes will end up growing much larger and 
radio luminosity falls much more quickly in the evolution of the 
source. Volume filling factors are increased from case [A] to ~ 0.4. 
The contribution to the unresolved CXB is 5 per cent. It may be 
the case that the typical density profile around FR II sources falls 
more steeply at higher radii, such as if the density profile is simi- 
lar to that of relaxed clusters dVikhlinin et aiT 2006). which follow 
the Navarro-Frenk-White halo profile (less ste ep than r~ 2 near the 
centre , more steep than r~ 2 at large distances) dNavarro et al J 19961 
Since most sources at redshifts z ^ 2 tend to fall below 
the radio flux limit early in their evolution before they reach Mpc 
scales (their fractional duration of visi bility compared to t he jet 
life time may be less less than 0.01-0.1) dBlundell et al J 19991) . and 
the total density of double-lobed sources depends on how quickly 
the sources fall below the flux limit, the density profile nearer the 
source rather than at large distances plays a more important role in 
determining the total population of double-lobed sources. 

If most sources were born over a broader span of time, which 
we investigate by setting Az = 1, the volume filling factor dur- 
ing the quasar era is slightly larger than estimated with our original 
parameter. This is a result of the fraction of sources that would 
be observable above the radio flux limit decreasing due to a wider 
spread in ages of sources at each redshift. With Az — 1, we find a 
volume filling factor of 0.07 at redshifts 2 and 3, assuming all the 
parameters are the same as in our original case [A], while the num- 
ber density and unresolved CXB contribution are comparable to the 
population described by the original birth function with Az = 0.6. 

Steepening the power-law index of jet energies from —2.6 to 
—3 has the effect of a larger number of sources being undetected 
in the survey and flattening it has the opposite effect. Compared 
to the original case [A] which predicts ~ 29 lobed-galaxies per 
square degree visible in the X-ray, of which ~ 4 are IC ghosts, 
above the X-ray flux limit of the CDFN survey, and predicts the 
radio sources contribute about 2 per cent of unresolved CXB and 
the volume filling factor of the lobes is ~ 0.03-0.02 at z — 2-3, if 
we increase the power-law index of jet energies from —2.6 to —3 
we find ~ 39 visible sources per square degree of which ~ 4 are 
IC ghosts, a 6 per cent contribution to the unresolved CXB, and 
a volume filling factor of the lobes is ~ 0.07-0.04 at z = 2-3. 
And if we decrease the power-law index of the jet energies to —2, 
we find ~ 16 visible sources per square degree of which ~ 2 are 
IC ghosts, a 1 per cent contribution to the unresolved CXB, and a 
volume filling factor of the lobes is ~ 0.01-0.004at z = 2-3. 



8 DISCUSSION AND CONCLUSION 

We have developed an analytic model to study the evolution of the 
radio luminosity and X-ray luminosity (due to IC scattering of the 
CMB) of FR II radio galaxies in order to quantify the abundance of 
actual and observable powerful double-lobed radio sources and IC 
ghosts (sources with jets turned off that still radiate in the X-ray). 
For a set of model parameters inferred from observations of radio 
sources (case [A]), which had tj = 10 8 yr, 7 m i n = 1, p = 2.14, 
po = 1.67 x 10~ 23 kg m~ 3 and /3 = 1.5, we predict ~ 29 lobed- 
galaxies per square degree visible in the X-ray, of which ~ 4 are 
IC ghosts, above the X-ray flux limit of the CDFN survey. The 
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CDFN survey found 167Ig 7 extended X-ray sources, which con- 
sist of both X-ray clusters and lobes. Thus a considerable fraction 
of these objects may turn out to be in fact powerful radio lobes 
based on our estimate. It is important to note that the region ob- 
served with the CDFN survey may not be a typical area of the X- 
ray sky because the survey was chosen to follow up on the Hubble 
Deep Field North observations, which looked at an area of the sky 
free from bright emissions in the visible, radio, infrared, ultraviolet 
and X-ray. It is also expected by our model that the powerful radio 
sources will contribute about 2 per cent of unresolved CXB, which 
is due to IC scattering of the CMB by the lobes and nuclei of in- 
dividual galaxies as well as clusters. The volume filling factor of 
the lobes is highest during the quasar era, at £ ~ 0.03 for 2 = 2. 
However, we also find that our predictions for volume filling factor 
are sensitive to the model parameters we use. Especially impor- 
tant could be the minimum injection energy, determined by 7 m j n , 
whose typical value is not well known. In addition, a higher value 
of 7min also improves the visibility of IC ghosts and can increase 
the observable ratio of IC ghosts to active sources by a factor of 3. 

For case [A] the volume filling factor is < 0.03, highest at 
2 ~ 2. The volume filling factor at the quasar era is about an order 
of magnitude above that of the present era. The volume filling factor 
declines only slightly from z — 2 to z — 3. The volume filling 
factor of dead lobes is much larger than that of active ones by at 
least an order of magnitude at all redshifts, but especially at low 
redshifts z < 1.5. 

The derived XLF at 1 keV of double-lobed sources (both 
active and ghosts) shows that they are more abundant at higher 
redshifts (2 > 1) than at lower redshift and are also generally 
more luminous in the X-ray at these higher redshifts (the peak 
of the XLF shifts to the right as redshift is increased). Addition- 
ally, at redshifts 2 > 2 the space densities of double-lobed radio 
galaxies become comparable to that of X-ray clusters at high lu- 
minosities (Z/ x ^ few x 10 43 erg s _1 ) - see Figure [Tol This is 
even true of just the high redshift IC ghosts at high luminosities 
(L x > few x 10 44 ergs -1 ). 

In this paper we have only considered FR II sources and 
have neglected FR I sources, which are much more n umerous 
in the local universe than FR lis. IWilman et ail d2008h creates 
a semi-empirical simulation of the extragalactic radio continuum 
sky out to redshift 2 = 20, and down to flux density lim- 
its of 10 njy, which i ncludes FR I and FR II simulated sources 
drawn from the RLF o f lWillottet all [2001) and shows that across 
all redshifts, FR I sources can be 5 orders of ma gnitude more 
abund ant that FR II sources. However, the RLF of Will ott et al.l 
d200 lh shows that FR I and FR II sources have comparable num- 
ber density during the redshift range (2 = 2 — 3), the quasar 
era, on which we are focus ing. Accessing simulation data prod- 
ucts of lWilmanetai] (2008) through the data base web interface 
(http://s-cubed.physics.ox.ac.ukl also agrees with there being a 
comparable number density during the quasar era. The relative con- 
tribution from FR I sources during the quasar era, for example to 
the volume filling factor, is modest. The volume of a typical FR 
I will be much smaller than the volume of an FR II source, since 
FR I seem to have a much higher efficiency i n converting je t thrus t 
to radio flux and are also much dimmer (see lBarai & Wiitj i2007t) 
and references therein). In assessing the importance of FR I ob- 
jects to the CXB, the question at hand is the integrated number of 
7 ~ 10 3 particles that the FR I produces during its lifetime. The 
number of particles ejected by an FR I will be less than that of an 
FR II object since the FR II sources have higher luminosity and jet 
power, and lower efficiency in converting beam power to radio flux. 



In addition, FR Is may have higher typical B-fields, which would 
reduce the number of particles in the lobe needed to produce the 
object's observed luminosity. The environ ments of an FR I object 
may be more magnetis ed. See, for example Laing et al.|j2008bl) and 
lLaing & Bridle! ([2008) where the FR I object 3C 3 l's magnetoionic 
environment is thought to have a B-field of 0.21 nT. This value is 
a lower limit to the B-field that 3C3 l's jet plasma is bathed in. The 
actual B-field in the lobes is likely higher since FR l's are typi- 
cally very lossy and local shear effects will only heighten the field 
strengths. While there are individual examples of FR Is that are 
exquisite ly well-modelled in t erms of velo c ity pro file, etc. (see for 
example, lLaing et alj J2008al) . lLaing et alj (2008b)) the modelling 
of typical FR Is to date is far behind that of FR lis. In a future 
work, it would be worthwhile to develop a model of FR I sources 
with comparable physical detail as the models for FR II sources to 
assess their contribution to the unresolved CXB. 

It is possible that typical double-lobed sources have mul- 
tiple outbursts of jet activity, while our model only accounted 
for a single outburst. I f sour ces show intermittent activity as in 
Reynolds & Begelmanl dl997h . active sources having fewer past 
outbursts will be brighter than those with a long history of episodes 
of jet ejection and hence more easily detectable. If during the jet 
lifetimes we assumed the source switches off a number of times 
then the lobes grow smaller in volume. Since sources tend to fall 
below the radio flux limit during the quasar era early in their life- 
times, the predicted total density of double radio sources will not be 
altered significantly by the inclusion of multiple outbursts (unless 
periods of jet activity are shorter than time at which the source falls 
below the radio flux limit). Therefore, the volume filling factor will 
be somewhat reduced. If, on the other hand, sources have recurring 
extended periods of jet activity of the order of the jet lifetimes we 
assumed, then lobes will grow larger still and the volume filling 
factor would increase. 

Our predicted total comoving density of active and non-active 
double radio sources at 2 = 1 and 2 is similar to the comoving 
density of galaxies above a lumi nosity Mk = — 25 and —24.5 in 
the galaxy luminosity function of Icirasuolo et al.1 J2010h . The pre- 
dicted total comoving density of FR II sources suggests that most 
massive galaxies in the universe may have an FR II type outburst 
at some point in their evolution. Using the correlation between K- 
band bulge luminosity and black hole mass for galaxies given by 
Marconi & Hun t (2003), the predicted total comoving density of 
active and non-active double radio sources is similar to the comov- 
ing density of all the brightest galaxies with su permassive black 
holes of masses > 10 87 M and 10 8 ' 4 Mrr jMcLure & Dunlod 
(2002) argue that the true radio loud AGN have black hole masses 
^ 10 8 ' 5 Mq. All these galaxies may at some point in their history 
have contributed to returning material into the IGM in large-scale 
outbursts. 

We have studied the actual and observable distribution of ac- 
tive double-lobed galaxies, IC ghosts and dead lobes. IC ghosts 
may comprise 10-30 per cent of double-lobed X-ray sources and 
could outnumber X-ray clusters at high luminosities. 

The volume filling factor of all lobes calculated from the wide 
range of model parameters investigated suggests that it is likely be- 
tween 0.001-1 during the quasar era, and the contribution to the 
unresolved CXB is 1-60 per cent. It is important to know the pa- 
rameters of a typical radio source well to be able to constrain the 
inferred distribution, as we see that our predictions have some ex- 
treme sensitivities to underlying physical parameters (§0. Having 
a significant fraction of the relevant volume of the universe per- 
meated by overpressured lobes of galaxies can trigger large scale 
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star for mation and seed filaments with magnetic fields, as summa- 
rized bv lGopal-Krishna & WiitJ < l200ll) . The volume filling factor 
at z = 3 is comparable to the volume filling factor at z = 2, 
meaning that radio lobes may have played an important role in 
star formation from an early time for a long-range period. Parame- 
ters of individual sources, such as 7 m i n , may significantly influence 
the cosmological predictions for volume filling factor and radiation 
contribution to the unresolved CXB, and therefore it is important 
to obtain better constraints on these parameters from observations. 
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